Magnetic properties of the low-dimensional spin-1/2 magnet a-Cu2As207 
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In this work we study the interplay between the crystal structure and magnetism of the py- 
roarsenate Qf-Cu2As207 by means of magnetization, heat capacity, electron spin resonance and 
nuclear magnetic resonance measurements as well as density functional theory (DFT) calculations 
and quantum Monte Carlo (QMC) simulations. The data reveal that the magnetic Cu-0 chains 
in the crystal structure represent a realization of a quasi-one dimensional (ID) coupled alternating 
spin-1/2 Heisenberg chain model with relevant pathways through non- magnetic As04 tetrahedra. 
Owing to residual 3D interactions antiferromagnetic long range ordering at Tn — lOK takes place. 
Application of external magnetic field B along the magnetically easy axis induces the transition to a 
spin-flop phase at Bsf ~1.7T (2 K). The experimental data suggest that substantial quantum spin 
fluctuations take place at low magnetic fields in the ordered state. DFT calculations confirm the 
quasi-one-dimensional nature of the spin lattice, with the leading coupling Ji within the structural 
dimers. QMC fits to the magnetic susceptibility evaluate Ji = 164 K, the weaker intrachain coupling 
Ji/Ji = 0.55, and the effective interchain coupling Jici/Ji = 0.20. 

PACS numbers: 71.20.Ps, 75.10.Jm, 75.30.Et, 75.50.Ee, 76.30.-v, 76.60.-k 



I. INTRODUCTION 



Crystals of composition X2Y20-J tend to adopt dif- 
ferent structures in which the Y atom may be octa- 
hedrally or tetrahedrally coordinated^. These systems 
have been found to exhibit interesting magnetic prop- 
erties ranging from a geometry induced spin frustration 
phenomena when X and Y are magnetic metal ions (see, 
e.g., Refi^), to a peculiar interplay between structure and 
magnetism in those compounds where Y is not mag- 
netic (see, e.g., Refsi^ii). Co, Ni, and Cu diphosphates 
and divanadates are interesting examples of those latter 
compounds. In their crystal structure two YO/^ tetra- 
hedra sharing an oxygen atom usually undergo a transi- 
tion from a staggered (180°) to a bent Y -0-Y config- 
uration (from thortveitite to dichromate type structure) 
with the decrease of temperature. A recent compara- 
tive band-structure study of copper (II) based pyrocom- 
pounds Cu2 1^207 ( Y=P, As, V)^ shows that despite hav- 
ing the same basic structural motif these materials rep- 
resent completely different magnetic models. Here the 
change of the non-magnetic central atom Y in the tetra- 
hedral anion groups FO4, e.g. from F = P to V, turns 
out to be crucial in determining the effective interchain 
exchange pathways depending on the electronic state {p- 
or d-type) of the central Y cation. Moreover, signifi- 
cant variations in the intrachain exchange coupling evi- 
dence the subtle effect of the Cu-O-Cu bond angle albeit 
the basic structural dimerization of CU2O6 plaquettes is 



present in all these isostructural compounds. Thus, ac- 
cording to theoretical predictions and experimental data 
available so far, one can expect that further change of 
the Y cation, might yield a completely new magnetic 
behavior in this series of materials. The validation of 
this conjecture would be of a specific importance for the 
development and verification of new magnetic models 
of low-dimensional Heisenberg quantum magnets which 
since decades stand in the focus of intense experimental 
and theoretical research in condensed matter physics^. 

Unlike pyrovanadates and even pyrophosphates, the 
experimental magnetic study of transition metal py- 
roarsenates, X2AS2O7 (X = Cu, Ni. Co, Mn) is scarce. 
Only for X2AS2O7 with X = Ni, Co, Mn, early bulk mag- 
netic susceptibility measurements and low-temperature 
neutron powder diffraction were reported by Buckley et 
al.^ whereas there seem to be no experimental results on 
the copper (Il)-based CU2AS2O7 pyroarsenate. There- 
fore, in view of the expected sensitivity of the magnetic 
behavior to the type of the Y cation we have carried out 
a detailed study of magnetic properties of single crystals 
of Q!-Cu2As207 with different experimental techniques 
such as magnetization, heat capacity, electron spin res- 
onance (ESR) and nuclear magnetic resonance (NMR) 
measurements, combined with density functional theory 
(DFT) calculations and quantum Monte Carlo (QMC) 
simulations. Indeed, we find that in contrast to other 
related isostructural compounds q;-Cu2As207 represents 
an experimental realization of coupled spin-1/2 alternat- 
ing antiferromagentic (AFM) Heisenberg chains where 
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FIG. 1. (Color online) The crystal structure of Q-CU2AS2O7: 
(a) - edge-sharing CU2O6 plaquettes form structural dimer 
units along the 6-axis. The dimers are linked by ASO4 anion 
groups forming a plane of coupled dimer chains; (b) - side view 
of the chain planes linked with interjacent AS2O7 anions. 

the ASO4 anion groups play an important role for the in- 
terchain magnetic exchange. The quasi-one dimensional 
magnetic behavior with short-range AFM correlations 
is dominating over a wide temperature range. At low 
temperatures, however, secondary magnetic exchange in- 
teractions become relevant, too, and three-dimensional 
long-range magnetic order evolves below Tn — 10 K. In 
the magnetically ordered state, the application of mod- 
erate magnetic fields i? ~ 1.7 T parallel to the crystallo- 
graphic 6-axis yields a field induced reorientation of the 
spin lattice. Further increase of the magnetic field pro- 
gressively suppresses spin fluctuations which are evident 
in the thermodynamic and magnetic resonance response 
at small fields. 



II. CRYSTAL STRUCTURE 

Synthesis and characterization of q;-Cu2As207 single 
crystals have been reported in Ref4. The a-Cu2As2 07 
modification crystallizes isotypically with a-Cu2P207 
and /3-CU2V2O7 structures, belonging to the monoclinic 
space group C2/c with a = 7.237(3), b = 8.2557(17), 
c = 9.780(3) A, /3= 111.03(2)°. The major structural 
feature of q;-Cu2As207 is the occurrence of non-planar 
CusOg plaquettes (Cu-O-Cu bond of 101.68(8)° and Cu- 
O bond lengths of 1.92-2.00 A) forming structural dimers 
along the crystallographic b direction^ (Fig- QJ^)- The 
ASO4 tetrahedra link these Cu-O-Cu dimers in a chain 
along the b-axis and also interconnect the neighboring 
chains in the chain plane (Fig. [1^). The chain planes 
are linked by interjacent pyroarsenate anions AS2O7 in a 
bent As-O-As configuration of 145.9(2)°. The 4-fold oxy- 
gen in-plane coordination of Cu is completed by the fifth 
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FIG. 2. (Color online) Magnetic susceptibility x ~ M/B of 
a a-Cu2As207 single crystal measured at _B = 1 T applied 
along the crystallographic b- and c-axis, respectively. 

oxygen atom from the neighboring plane with a much 
larger bond length of 2.27 A which therefore should only 
weakly contribute to magnetism. 

III. RESULTS AND DISCUSSION 
A. Magnetization 

Magnetization measurements on a single crystal of a- 
CU2AS2O7 were carried out in a VSM SQUID magne- 
tometer (Quantum Design) applying the magnetic fields 
parallel (_B|| 6-axis) and perpendicular (i3|| c-axis) to the 
propagation direction of the CU2O6 plaquettes, respec- 
tively. Fig. [5] shows the static magnetic susceptibility 
X = M/B in the temperature range 2 - 350 K at B = 1 T. 
For both field directions, the susceptibility data exhibit a 
characteristic broad maximum at about T ~ 105 K which 
is typically associated with ID short-range magnetic cor- 
relations due to predominant intra-chain magnetic inter- 
actions. This value suggests thus a rough estimate of the 
expected energy scale of the leading exchange constants. 
At a much smaller temperature Tn — 10 K, there is an 
upturn in x(T)||c and a concomitant kink in x(T)||b evi- 
dencing the onset of long-range AFM spin order (see also 
figures [5] and [TU)) . Such a feature and the rapid decrease 
of x(T)||& down to T = 2 K imply the crystallographic 
6-axis being the magnetically easy axis. 

Further x(T)||^ data at various constant applied mag- 
netic fields are shown in Fig. |3l At low field _B < 2 T. the 
magnetic specific heat which is proportional to 9(xT) / dT 
exhibits a pronounced A-like anomaly (Fisher's heat 
capacity)^iiS corroborating the scenario of long range 
magnetic ordering at Tn = 10 K (Fig. [3l inset). In- 
terestingly, the kink in x{T) at Tn significantly changes 
its shape upon application of magnetic field _B > 2 T. 
The decrease of x{T) below Tn turns into a steep upturn 
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FIG. 3. (Color online) Temperature dependence of the mag- 
netic susceptibility x(.T)- Inset: The pronounced anomaly in 
the Fisher's heat capacity d{xT)/dT at B = IT along the 
b-axis implies long-range magnetic order below Tn — 10 K. 

at higher fields (Fig. [3]). This transformation should be 
indicative of an additional magnetic phase transition at 
a critical field in between these values and the develop- 
ment of a field-induced magnetic phase. Such a behavior 
is indeed clearly visible in the magnetic field dependence 
of the magnetization M vs. B\\b (Fig. |4]). 

Field dependent measurements of magnetization for 
i3||6-axis in Fig. 2] give evidence of a first-order spin- 
fiop transition (SF) below Tn at Bsf ^ 1-7 T (2 K), in 
agreement with the x(T) data. The spin-flop arises due 
to the gain in magnetic anisotropy energy of the system 
above Bsf considering the AF collinear arrangement 
of Cu^+ spins parallel to the b-axis at lower fields. The 
magnetization AI{B) data for i3||c-axis in the field range 
up to 7 T (not shown) are consistent with a magnetically 
hard c-axis, where no anomaly is observed below 12 K. 



B. Heat capacity 

Fig. [5] shows the molar specific heat at B = and in 
different external magnetic fields applied parallel to the 
&-axis. At i? = 0, the as measured data do not exhibit 
a pronounced anomaly at Tn despite the marked peak in 
the d{xT)/dT (Fig. [3]) clearly indicates anomalous mag- 
netic entropy changes associated with the onset of long 
range AFM order. We conclude that, in zero magnetic 
field, the magnetic entropy changes at Tn are small com- 
pared to the phononic ones. In order to highlight this, 
we have subtracted an arbitrary polynomial background 
function f(T) = aT + bT'^ from the experinrental data 
(cf., e.g.^). The resulting data Acp are shown in the 
bottom panel of Fig. [5] At a temperature very close 
to Tn derived from the susceptibility measurements the 
data indicate a kink in Acp, at i? = 0, which is associated 
with a specific heat anomaly Acp < 20 mJmol~^K~-'^. 
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FIG. 4. (Color online) Magnetic field B dependence of the 
magnetization M for B || 6-axis at different temperatures 
slightly above and below Tn. CP and SF denote the pre- 
sumably collinear phase in the low-field regime and the spin- 
flop phase at high magnetic fields. The spin-fiop transition 
at _BsF — 1.7 T (T = 2K) is associated with the magnetic 
anisotropy in the system. 
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FIG. 5. (Color online) Upper panel: Temperature depen- 
dence of the specific heat capacity Cp at B = and 9 T; lower 
panel: Experimental data in different magnetic fields -B||&- 
axis after subtraction of an arbitrary polynomial background 
function /(T) = aT+bT'^. indicates the temperature of 
the anomaly in d{xT)/dT at small fields (cf. Fig. |3j. 
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This kink is hardly affected by small magnetic fields. To 
be specific, at B = 5T the kink is fomid at T^^^(5 
T)= 10.2 K. However, at magnetic fields B > 5T the 
peak increases in amplitude and shifts to higher temper- 
atures. The growth of the peak suggests the increase 
of entropy changes associated with the phase transition, 
i.e. the evolution of a jump-like or a A-like specific heat 
anomaly. At i? = 9 T, the peak height of the specific heat 
anomaly amounts to ~ 115 mJmol~^K~^ and the peak 
maximum occurs at T^(9 T)= 11.9 K. The experimental 
data hence prove a field induced increase of the mag- 
netic entropy changes associated with the onset of long- 
range magnetic order. The vanishing size of the specific 
heat anomaly in small magnetic fields, i.e. the absence 
of significant entropy changes associated with the onset 
of AFM order, either implies only weak spin order below 
Tn or the presence of strong spin correlations well above 
the long-range ordering temperature. The fact that ap- 
plication of moderate external magnetic fields instead of 
suppression of AFM order yields larger entropy changes 
associated with the phase transition corroborates the lat- 
ter scenario. Hence, the specific heat data suggest strong 
AFM spin fiuctuations at low fields. We note that from 
the present data it is difficult to conclude if additional 
small shoulder-like features of the Acp(r, B = 9T) curve 
at ~ lOK and ^ 13.3 K are of intrinsic nature or arise 
due to the imperfect subtraction of the nonmagnetic con- 
tributions from the total specific heat. 

A tentative phase diagram which has been constructed 
from the magnetization and specific heat data is pre- 
sented in Fig.ini In particular, it shows that the spin-flop 
field which signals the transition from the AFM collinear 
phase (CP) to the spin-flop phase (SP) is nearly in- 
dependent on temperature. For a discontinuous phase 
transition, the slope of the phase boundary is associ- 
ated with the ratio of the magnetization changes AM 
and the entropy changes AS at the phase transition 
according to dT/dB = ~AA!/AS (Clausius-Clapeyron 
relationi^). Hence, while AM > our observation of 
dB/dT ~ implies the absence of significant entropy 
changes at the spin-flop transition, as expected for a spin- 
reorientation in the AFM ordered phase. Quantitatively, 
at T = 1.7 K, we extract dT/dB « 16.2 K/T and AM = 
2.4- 10-3 /^B/f.u. which yields ASsf = 0.83 mJ/(Mol-K). 
Consquently, the degree of spin order at high fields/low 
temperature in the spin-flop phase is very similar to that 
at zero field/low temperature. Applying the above re- 
lation at T = 5 K and T = 7K yields decreasing en- 
tropy changes at the spin flop transition AS'sf(5K) = 
0.64mJ/(Mol-K) and AS'sf(7K) 0.45mJ/(Mol-K) as 
already qualitatively visible from the rather constant 
slope of the phase boundary while AM decreases with 
temperature. 

In contrast, the phase boundary between the param- 
agnetic and the AFM CP/SF phases implies dT/dB « 
for S < 5 T and dT/dB > for higher flelds. Our data 
do not allow to judge the exact nature of the speciflc 
heat anomaly and to quantitatively estimate the entropy 
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FIG. 6. (Color online) Magnetic phase diagram for a- 
CU2AS2O7 single crystals according to magnetization and spe- 
cific heat measurements when applying field along the mag- 
netically easy 6-axis. 



changes or the specific heat jump associated to the phase 
transition. Firstly, dT/dB « thermodynamically cor- 
responds to only tiny entropy changes at Tn(B < 5 T) if 
the Ehrenfest equationi^ (related to a continuous phase 
transition) and the moderate kink in dM /dT at the tran- 
sition are considered. This indeed agrees to our specific 
data. However, the pronounced positive slope at the PP- 
SF transition appears in qualitative contrast to what is 
expected for an ordinary AFM phase transition. 



C. Electron Spin Resonance (ESR) 

ESR at a frequency of ~ 10 GHz has been measured 
with a commercial X-band spectrometer (Bruker KMX). 
High-field high-frequency ESR (HF-ESR) measurements 
have been carried out with the home-made spectrome- 
ter on the basis of the Millimeterwave Network Analyzer 
from AB Millimetre and a 15 T superconducting magnet 
from Oxford Instruments^^ A well defined ESR signal 
of a Lorentzian shape associated with the resonance of 
Cu^+ {3d^, S ~ 1/2) ions was observed in a broad range 
of excitation frequencies and temperatures. 

The dependence of the position of the ESR line B^^s 
on the direction of magnetic field within the crystallo- 
graphic be plane has been studied at ~10 GHz at room 
temperature (not shown). The g-factor calculated from 
this dependence as 5 = hv / n^B^^s reveals an anisotropy 
typical for Cu^"*" ions in a square planar coordination: ^^1^''' 
g\\ ~ 2.07 {B\\b - along the chains) and g±_ ~ 2.18 {B\\c). 
Here h is the Planck constant and /ie is the Bohr mag- 
neton. 

Temperature dependent HF-ESR measurements were 
carried out aX v = 340 GHz. For magnetic fields applied 
along the magnetic easy axis (_B||6), the ^-factor is almost 
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temperature independent above ~60 K (Fig. [T]) close to 
a value of (?|| ~ 2.09 which corresponds to B^cs — 11.6 T. 
A slight discrepancy to the 10 GHz data (gy ~ 2.07) is 
likely due to a slight misorientation of the sample in the 
HF-ESR probe head. At lower temperature the g-factor 
drops with decreasing temperature indicating the onset 
of short-range magnetic ordering when approaching Tjy. 
The shift of the resonance field in this region concomi- 
tant with the broadening of the line is apparently due to 
the growth of local quasi-static internal fields and the in- 
crease of the spin correlation length as precursors of the 
long range magnetic order. 

Frequency dependent measurements at T = 20 K > Tn 
show a linear scaling between v and B (Fig. [8]) yielding a 
g-factor g ~ 2.07 from the slope dv/dB. Extrapolation 
of the data to zero field reveals no energy gap for the 
resonant excitation as expected for an S* = 1 /2 system in 
the paramagnetic state. An opening of a gap at B = is 
however expected in the magnetically ordered state due 
to the magnetic anisotropy. For a simple two-sublattice 
collinear antiferromagnet the anisotropy gap A(_B = 0) 
is related to the spin-fiop field i?sf as A = BsigjiB/hi^ 
With BsF =i 1.7 T one obtains A = 50 GHz. By apply- 
ing the field along the easy 6-axis two resonance modes 
should develop as sketched in Fig. |S1 which collapse at 
BgpJ^. At i? > BsF one resonance branch is expected 
which approaches the paramagnetic ^{B) line in strong 
fields, as shown in Fig. [51 Indeed at high fields a well de- 
fined resonance mode has been detected at T = 4 K < Tn 
which follows closely the predicted field dependence (Fig. 
[5]) . Surprisingly it fades by approaching the field of 5 T 
and no resonance signals could be found at smaller fields. 
Such peculiar absence of resonance modes in the low field 
regime is in a remarkable correspondence with the van- 
ishing of the specific heat anomaly below 5 T which 
again points at the pronounced spin fluctuations in the 
magnetically ordered state which can be suppressed by 
magnetic field. 



2.10 



2.05 



o 2.00 
> 



1.95 



1.90 




1.2 
1.0 
0.8 

r; 
0.6 I 
? 

0.43 
0.2 
0.0 



40 80 120 160 200 240 
T(K) 

FIG. 7. (Color online) Temperature dependence of the ESR 
parameters with B\\b at 340 GHz. Main panel: effective g- 
factor and linewidth, upper and lower curve respectively. In- 
set: The shift of the ESR resonance line at low temperature 
reflects the increasing internal magnetic field in the vicinity 
of the ordered state. 
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D. Nuclear Magnetic Resonance (NMR) 



FIG. 8. (Color online) Frequency versus resonance field di- 
agram of the ESR signal in the ordered and paramagnetic 
states for -B || 6. At 4K a single resonance mode is observed 
above ~5 T (solid circles). Solid thin line corresponds to the 
paramagnetic resonance branch v — gyL^Brca/h. Solid thick 
lines show the expected AFM resonance branches for uniaxial 
anisotropy with B \\ easy b-axis.— 



NMR on ""^As nuclei was measured with a Tecmag 
pulse solid-state NMR spectrometer with a 9.2 T super- 
conducting magnet from Magnex Scientific. The NMR 
spectra were obtained by measuring the intensity of the 
Hahn echo versus magnetic field. The Ti relaxation time 
was measured with the method of stimulated echo. 

Representative ^^As NMR spectra for the easy axis h- 
direction of the magnetic field are shown in Fig. [HI The 
temperature dependence of the ^^As spin lattice relax- 
ation rate Tj~^ divided over T is shown in Fig. [TUl On 
this figure it is compared with the T-dependence of the 
static susceptibility x{T). The fact that both (TiT)"! 
and x{T) follow closely the same temperature behavior 
suggests an appreciable coupling of the ^^As nuclei to the 



Gu spin system. In this case the nuclear spin lattice re- 
laxation is usually governed by magnetic fiuctuations in 
the electronic spin system and (TiT)~^ is proportional 
to the imaginary part of the dynamic susceptibility x": 

(rrT)-V^Ai(g».^^ (1) 

9 

Here A is the hyperfine constant, q is the wave vector and 
w is the NMR frequency. Note a sharp peak of (TiT)~^ 
at the ordering temperature Tn manifesting a critical be- 
havior of Gu spin fluctuations (Fig. (TUl) . At T > Tn 
the ^5 As NMR spectr um consists of three quadrupole- 
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FIG. 9. (Color online) '^^As NMR spectra (dots) measured at 
30 K and 4 K for the magnetic field direction parallel to the 
6-axis. Solid lines represent simulated spectra (see the text). 



on the basis of Hamiltonian ^ with the same values of 
vq and 77 assuming two magnetically nonequivalent As 
sites (Fig. [5]). The modelling yields the magnitude of 
the internal field at the ^^As amounting to Bint = 0.3 T 
which is directed approximately perpendicular to the b- 
axis. The substantially large value of Bint suggests a 
transferred nature of the internal field due to the overlap 
of electronic orbitals. 

The finding that Bint -L b enables a conclusion that 
though the NMR has been measured at fields ^ 7 T 
substantially stronger that the spin-flop field Bsf — 1.7 T 
the Cu spins in the chains are still oriented predominantly 
perpendicular to the easy axis. This suggests the AFM 
exchange interaction between the Cu spins to be much 
stronger compared to the Zeeman energy which renders 
the spin canting in a field of ~ 7 T very small. This 
complies with the small magnetization values shown in 
Fig.H 
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FIG. 10. (Color online) Temperature dependence of the ^^As 
spin lattice relaxation rate {T^^T)~^ in comparison with the 
static susceptibility. Thin solid line is the guide for the eye. 



split lines as expected for nuclei with the spin 7 = 3/2 
(Fig. [H]) . The spectrum has been modeled on the basis of 
the following Hamiltonian>ii 
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The first term here describes the quadrupole interaction 
and the second term the Zeeman interaction. A good 
agreement with the experiment has been achieved with 
the value of the quadrupole frequency j/q = 23 MHz and 
the asymmetry parameter = 0.1. In the magnetically 
ordered state (T < Tn) the spectral lines experience a 
splitting into two components suggesting the occurrence 
of two nonequivalent magnetic environments of the ^^As 
nuclei. This can be related to the occurrence of two spin 
sublattices in the ordered Cu chains which exert alter- 
nating in sign internal field at the As sites. The low 
temperature NMR spectra can be successfully modeled 



E. Density functional band structure calculations 

DFT calculations have been performed using the full- 
potential code FPLO-8.50-32fi§. For the exchange and 
correlation potential, the scalar-relativistic parameteriza- 
tion of Perdew and Wang has been usedJ^ The fc-mesh 
of 1152 points (320 in the irreducible wedge) has been 
adopted for the nonmagnetic calculations. The transfer 
integrals ti have been evaluated based on a one-orbital 
tight-binding model, parameterized using Wannier func- 
tions. Strong electronic correlations are accounted for on 
a model level by resorting to a Hubbard model. At half- 
filling, the low-energy (magnetic) sector of this model 
is effectively described by a Heisenberg model.— Anti- 
ferromagnetic exchange integrals J^^^ have been deter- 
mined based on transfer integrals ti, adopting the second- 
order perturbation theory expression jf^^^ = Atf/Ucff- 
Alternatively, the local (spin) density approximation 
[L(S)DA]+C/ has been used to evaluate the values of to- 
tal exchange, comprising AFM as well as ferromagnetic 
(FM) contributions: Ji = J^^^ + Jf^. The fc-meshes 
used for LSDA+J7 supercell calculations number 12, 18 
or 32 points, depending on the supercell size. The results 
were accurately checked for convergence. 

Typical for cuprates, the valence band of a-Cu2As207 
is dominated by Cu and O states (Fig. [TT|) . The nonzero 
density at the Fermi level evidences a metallic ground 
state, in contrast to the olive color of a-Cu2As207 crys- 
tals. This discrepancy, typical for the LDA, originates 
from insufficient treatment of strong electronic correla- 
tions, intrinsic for the 3dP electronic configuration of 
Cu^^. Despite the inappropriate treatment of electron 
interaction terms, LDA correctly accounts for kinetic (ex- 
change) processes, and thus can be used to derive the 
relevant orbitals and couplings. This way, we find that 
the magnetic properties of q;-Cu2As2 07 are governed by 
a Cu 3d orbital of the x'^ — symmetryi^ For these 
states, the Wannier functions are constructed^^. 
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FIG. 11. (Color online) LDA total and atomic-resolved den- 
sity of states of a-Cu2As207. The main contribution to the 
valence band is due to Cu and O states (full orange and dashed 
red line, respectively). The Fermi level is at zero energy. 
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FIG. 12. (Color online) Microscopic spin model for a- 
CU2AS2O7. Magnetic layers are formed by alternating Ji — J'l 
chains running along b that are coupled by two inequivalent 
couplings Jici and Jics. The leading coupling Ji corresponds 
to the magnetic exchange within the structural CU2O6 dimers. 



The analysis based on Wannier functions yields four 
relevant (> 30 meV) transfer integrals. The leading terms 
are intra- and inter dimer couplings ti and t'l, forming al- 
ternating chains running along b. In addition, the chains 
are coupled together by two inequivalent transfer inte- 
grals iici and tics shaping a 2D spin lattice (Fig. [T^ . 
Numerical values of the relevant U terms are given in 
Table |T] (third column) . 

Mapping onto a Hubbard model and subsequently onto 
a Heisenberg model yields the values of j/^^^ (Table [l| . 
Considering Tmax = 105 K from our x(T) measurements, 
the large J^^^ = 298 K resulting from a typical value 
Ucs = 4.5 eV hints at a substantial FM contribution to 
this exchange. To evaluate the values of the total ex- 
change Ji, we perform LSDA+t/ calculations. This cal- 
culational scheme adds correlations on top of the L(S)DA 
in a mean-field way. The implementation of the method 
relies on several parameters, namely, the on-site intra- 
orbital term Ud, the on-site inter-orbital term J^, and the 
particular way of accounting for the part of correlation 
energy already present in LDA. The latter term is typi- 
cally referred as the double-counting correction (DCC), 
and two most widely used schemes are the around-mean- 



TABLE I. Leading couplings in a-Cu2As207; Cu-Cu dis- 
tance d (in A), transfer integrals ti (in meV) and AFM ex- 
change integrals jf^^^ (in K), calculated adopting [/off = 4.5 
eV. Total exchange integrals Ji (in K) are evaluated from 
LSDA+t/ calculations adopting AMF {Ud = 6.5 eV) and FLL 
[Ud = 9.0 eV) DCC. The paths are depicted in Fig. [H 



path 


d 


U 


rAFM 


Ji, AMF 


J^, FLL 


Ji 


3.086 


170 


298 


167 


188 


Ji 


5.239 


104 


112 


125 


102 


Jicl 


6.437 


76 


60 


44 


42 


Jic3 


6.465 


77 


61 


46 


42 



field (AMF), which uses a uniform orbital occupation as a 
reference, and the fully localized limit (FLL), referring to 
integer (0 or 1) occupation numbersi^ For [7^, we used 
the range 6.5±1 eV within AMF and 8.5±1 eV within 
FLL. keeping J^j = 1 eV constant. 

Representative values of the resulting J,; are given in 
Table HP^ In agreement with the Goodenough-Kanamori 
rules,^^ the Cu-O-Cu angle of 101.7° leads to sizable 
FM contribution, and consequently, to a substantial re- 
duction of Ji. However, in contrast to the isostruc- 
tural compounds Q!-Cu2P207 (Ref. Q) and /3-CU2V2O7 
(Ref.i), this FM contribution does not suffice to reduce 
Ji strongly, thus it remains the leadiirg coupling in a- 
CU2AS2O7. FM contributions to other relevant couplings 
are substantially smaller due to their long-range nature. 

The interchain coupling topology in a-Cu2As207 is 
similar to Q!-Cu2P207, but contrasts to /3-CU2V2O7. As 
discussed in Ref. , this essential difference originates 
from the nature of the atom residing in the center of an- 
ionic tetrahedra. in particular, whether it is a p-element 
(P, As) or a d element (V). The former case favors the 
Cu-O-O-Cu connections via single tetrahedra, i.e., J{ci 
and Jic3. 

To summarize the results of the band structure calcu- 
lations, the microscopic magnetic model of Q!-Cu2As207 
features alternating spin chains, with both nearest- 
neighbor couplings AFM, and the dominant coupling Ji 
within the structural CU2O6 dimers. The chains are cou- 
pled by two inequivalent interchain couplings, forming a 
2D nonfrustrated nrodel. 



F. QMC model studies 

DFT calculations typically provide numerically accu- 
rate spin models, with surprisingly small errors for the 
individual exchange integrals Ji. Aiming to evaluate 
a quantitative microscopic magnetic model with even 
higher precision, we perform QMC simulations of the ob- 
served quantities and subsequently compare them to the 
experiments We first attempt to describe the magnetic 
susceptibility of a-Cu2As207 for B\\b using the param- 
eterized solution32i for the alternating Heisenberg chain, 
valid in the temperature range O.lJifcg^ <T< iJik^^. 
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The fit exhibits good agreement with the experimental 
curve down to 40 K, yielding Ji = 174 K and J{ = 112 
K that conform to our LSDA+J7 results. However, the 
fitted value of factor amounts to only 1.94, strongly 
underestimating the intrinsic value of 2.07 derived from 
ESR. 



Next, we take the interchain couplings into account. 
This way, we calculated the reduced magnetic suscepti- 
bility X* using QMC simulations of the spin Hamilto- 
nian, by considering finite lattices of = 2048 sites (32 
coupled chains of 64 spins each) with periodic boundary 
conditions. We used the stochastic series expansion algo- 
rithm implemented in the code LOOpSS from the software 
package ALPSi^ Low statistical errors (below 0.1 %) are 
ensured by using 50 000 sweeps for thermalization and 
500 000 steps after thermalization. The resulting simu- 
lated x*{T*) dependencies were fitted to the experimen- 
tal x(T) curve using the expression: 



X{T) ^ I, J -X [t—TI+^o, (3) 



where xo is the temperature-independent contribution, 
and Na and /is are the Avogadro constant and the Bohr 
magneton, respectively. 



As a starting point, we adopt the ratios Ji : J( : Jici = 
1 : 0.75 : 0.25 from the LSDA+J7 AMF calculations (Ta- 
ble U last-but -one column). Although the resulting fit 
reasonably agrees with the experimental curve, small de- 
viations evidence notable differences between our micro- 
scopic model and the experimental data. To adjust the 
model parameters, we varied J(/Ji and Jici/Ji ratios 
in a broad range, while the two relevant interchain cou- 
plings Jici and Jic3 were kept equal. This way, the best 
fit is attained for Ji : J{ : Jid = 1 : 0.55 : 0.20, in al- 
most perfect agreement with the ratios 1 : 0.54 : 0.22 
from the FLL DCC of LSDA+C/ (Table H last column), 
with the absolute value of Ji = 164 K only 13 % smaller 
than its LSDA+C/ FLL estimate. Even more remarkable, 
the QMC fit yields g=2.06, in excellent agreement with 
the experimental ESR estimate of 2.07 for B\\b. There- 
fore, the QMC simulations reveal crucial importance of 
interchain couplings for a proper description of the ther- 
modynamic properties of Q!-Cu2As207. 




100 200 



r(K) 

FIG. 13. (Color online) QMC fits to the experimental (exp.) 
magnetic susceptibility. The Ji : J[ : Jid ratios are given in 
the legend. The difference plots (inset) indicate the relative 
goodness of fit (Ax/x)^ for the two parameter sets of Ji : J[ : 

Jici- 



To investigate the role of interchain couplings for the 
magnetically ordered ground state, we calculate diago- 
nal spin correlations {SqS§^) for < R <6 along different 
paths on the spin lattice (Fig. [T31 right panel). First, 
we address the correlations between the alternating Ji- 
J'l chains (Fig. [TH dashed and thin solid lines). As ex- 
pected, Jici and Jic3 give rise to substantial interchain 
correlations, and thus favor a quasi-2D magnetic model. 
However, the correlations along the Ji-J( alternating 
chains (Fig. [T31 thick line) nearly coincide with the spin 
correlations of an alternating Heisenberg chain featuring 
the same ratio (0.55) of the nearest-neighbor couplings 
(Fig. [131 dotted line). Therefore, the interchain cou- 
plings, despite their crucial role for long-range magnetic 
ordering, have only a minor impact on the key element 
of the magnetic model — Ji^J'i alternating spin chains. 
As a result, even the purely ID Ji-J[ model reasonably 
fits the experimental magnetic susceptibility, although 
the underestimated g value and the magnetic transition 
at Tn evidence notable deviations from such a simplified 
model. 
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FIG. 14. (Color online) Left panel: diagonal spin correla- 
tions {5o5'r)(R) in the magnetic ground state of CU2AS2O7. 
The correlations along J\-J[ bonds (thick solid line, empty 
squares) closely resemble those of an alternating Heisenberg 
chain with J(:Ji=0.55 (dotted line). In contrast, the cor- 
relations between the Ji~J[ chains are substantially smaller 
(dashed and thin solid lines). Right panel: 0-R pathes in the 
spin lattice of CU2AS2O7 used for the {5'oSr^)(R) plot in the 
left panel. 



G. Structure-property relationships 

Since continuous symmetries can not be spontaneously 
broken at finite temperature neither in one nor in two 
dimensions)^ the observed long-range AFM order be- 
low Tn = 10 K can be accounted for only by extending 
the model with an effective coupling between the mag- 
netic layers. Considering the small size of this coupling, 
its accurate estimate based on DFT results is at best 
challenging. In addition, the propagation vector of the 
magnetically-ordered structure is not known, impeding a 
conclusive analysis of the interchain coupling regime on a 
microscopic level. Therefore, the magnetic ordering tem- 
perature as well as the propagation vector are beyond the 
scope of the present study. 

Magnetic layers formed by coupled structural dimers 
are a common element in spin-^ systems based on Cu^+ 
[a-Cu2P2 07 (Ref. i) and /3-CU2V2O7 (Ref. |i)]as well 
as [(VO)2P207 (Ref. [H), CSV2O5 (Ref. M), and 
VOSeOa (Ref. 134!)]. A key feature of their microscopic 
magnetic models is the minor role of the coupling Ji 
within structural dimers, compared to the stronger cou- 
pling running via the double bridge of nonmagnetic an- 
ions between the structural dimers {J[ in our present 
notation) The reduction in Ji is generally understood 
as the suppression of Cu-O-Cu or V-O-V superexchange 
due to a bond angle slightly exceeding 90°. In con- 
trast, Q!-Cu2As207 stands out against this trend, since 
the dominating coupling Ji is restored unexpectedly. The 
enhanced value of Ji ~ 164 K is only partially related to 
the change in the Cu-O-Cu bond angle that amounts to 
101.7°, compared to 100.4° in a-Cu2P207 (Ji ~ 30 K)^ 
and 98.7° in /3-CU2V2O7 (Ji < 12 K)A An efficient mech- 
anism affecting the coupling could be the polarization of 



particular oxygen 2p orbitals by As states, as proposed 
for 06"*+ in CuGeOsi^ A microscopic investigation of 
other Cu^+ arsenates will help to clarify this issue. 

The coupling regime of q;-Cu2As207 underscores the 
wide diversity of spin lattices arising from the quasi- ID 
arrangement of structural dimers. Apart from the afore- 
mentioned Cu^"*" and V^"*" compounds, a similar struc- 
tural feature is found in Cu2(PO.^)2CIl2 (Ref. [37h and 
VOHPO4 • O.5H2O (Ref. [13. Altogether, these com- 
pounds exhibit remarkably different ground states, di- 
mensionality, and even frustration regimes: compare, for 
example, the spin gap in the quasi-lD (VO)2P207 to the 
long-range magnetic order in the quasi-2D /3-CU2V2O7, 
and further to the quasi-2D frustrated Cu2(P03)2CH2 
that again features a singlet ground state. q;-Cu2As207 
takes its own distinct position in this series as a quasi- 
lD nonfrustrated system exhibiting a magnetic ordering 
transition at the notably low T^/ Ji ~ 0.06. 



IV. CONCLUSIONS 

The magnetic properties of a-Cu2As207 single crys- 
tals have been investigated using several experimental 
methods (static susceptibility, magnetization, heat ca- 
pacity, ESR and NMR measurements) , as well as density- 
functional band-structure calculations and QMC simula- 
tions. The results show that q;-Cu2As207 is a spin-1/2 
low-dimensional system with the crystallographic 6-axis 
as the magnetically easy axis. The intrachain interac- 
tions are reflected in ID short range AFM correlations 
below ^100 K and short range magnetic ordering be- 
low '^60 K. A significant interchain coupling contribu- 
tion drives the system into an AFM long range order 
at Tjv ~ 10 K. We find that a second phase transition 
of spin-flop type takes place at T < Tjv for Bsf — 1-7 
T (2 K), most likely related to the gain in magnetic 
anisotropy energy. The coUinear AF magnetic ordering 
below Bsf confirmed by a well defined A-like anomaly 
in the magnetic specific heat quantity. On contrary, nei- 
ther specific heat shows the expected A-like specific heat 
anomaly below '^5,T, nor ESR signals are detectable in 
this low field range suggesting the occurrence of strong 
spin fluctuations at small fields. The occurrence of the 
AFM resonance mode and the splitting of the ^^As NMR 
signals observed in the high field regime suggests the ro- 
bust AFM spin order in the spin fiop phase above '-^5T. 
The DFT calculations yield a reliable determination of 
all relevant magnetic couplings. A combined analysis of 
experimental and theoretical data enables a conclusion 
that, in contrast to other isostructural compounds, the 
pyroarsenate q;-Cu2As207 can be effectively described by 
a coupled alternating Heisenberg chain model with the 
leading intrachain coupling placed within the structural 
dimers Cu206- The topology and the strength of ex- 
change couplings determined from DFT calculations ren- 
der a spin frustration scenario unlikely and rather put 
forward low-dimensional nature of the spin system. We 
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conclude that indeed the electronic state of the central 
ion in the non-magnetic ASO4 side groups plays a cru- 
cial role for determining the relevant interchain pathways 
which makes a-Cu2 AS2O7 distinct in its magnetic proper- 
ties from other representatives of this class of compounds. 
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